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Abstract. In this work we investigate the material point model (MP-model) 
and exploit the geometrical meaning of the "entropy form" introduced by 
B.Coleman and R.Owen ([?]). We analyze full and partial integrability (close- 
ness) condition of the entropy form for the model of thermoelastic point and 
for the the deformable ferroelectric crystal media point. We show that the ex- 
tended thermodynamical space introduced by R.Hermann and widely exploited 
by R. Mrugala with his collaborators and other researchers, extended possibly 
by time, with its canonical contact structure is an appropriate setting for the 
development of material point models in different physical situations. This 
allows us to formulate the model of a material point and the corresponding 
entropy form in terms similar to those of the homogeneous thermodynamics, 
1261 . Closeness condition of the entropy form is reformulated as the require- 
ment that the admissible processes curves belongs to the constitutive surface 
S of the model. Our principal result is the description of the constitutive 
surfaces of the material point model as the Legendre submanifolds S5 (equi- 
librium submanifolds of homogeneous thermodynamics) of the space V shifted 
by the flow of Reeb vector field. This shift is controlled, at the points of Le- 
gendre submanifold Eg by the entropy production function a. 
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1. Introduction. 

The goal of this work is to investigate the material point model (MP-model) and 
exploit the geometrical meaning of the "entropy form" introduced by B.Coleman 
and R.Owen ([7]) and, later on, apphed to different physical systems in plISllQllTOl 

In their work, B.Coleman and R.Owen developed the basic dynamical scheme of 
thermodynamical systems postulating the basic properties of states and processes 
in a thermodynamical systems, introduced the abstract notion of action, Clausius- 
Duhem inequality etc. 

Geometrical structure of their model was later on reformulated in terms of the 
bundle theory in [8|[9l[T0]. 

The second part of their work is devoted to illustrations and applications of their 
scheme to the theory of " simple material elements" , including the " elastic points" . 
The authors introduced the "entropy form", defined in the appropriate state space 
(extended by time variable) and determining the change of entropy produced by 
the flux along the process defined by the curve in the state space. Entropy form 
is constructed starting from the balance of entropy of continuum thermodynamics. 
Applying the energy balance and some kinematical relations for the rate of change 
of the involved dynamical variables, one rewrites the entropy increase due to the 
flux as the integral of a 1-form rj' in the state space. Entropy function is then defined 
as an upper potential of the entropy form. Such a potential exists, in particular, 
in the case where entropy 1-form is closed. Closure conditions provides some set 
of constitutional relations for participating fields. The study of the entropy form 
and of the corresponding constitutive relations given by the condition of closeness 
of the entropy form for different thermodynamical systems were done in the works 
[i [91 [TOl [12] and some others. 

In this work we would like to analyze the geometrical meaning of the entropy 
form and that of its integrability. 

In the first part we revisit the " basic model" of the thermoelastic material point 
of Coleman-Owen in order to determine which part of the constitutive relations 
for the dynamical system of MP-model can be obtained from different integrability 
conditions of the entropy form and what information should be added from the 
continuum thermodynamics in order to construct the closed dynamical system of 
MP-model. We also reformulate the "material point-entropy form" model both for 
the " basic model" of thermoelastic point and for the deformable ferroelectric crystal 
media 12J, in terms of extended thermodynamical phase space. We study integrabil- 
ity conditions of the entropy form in the way similar to the study of Coleman-Owen 
simple model. We determine the constitutive part of the integrability conditions 
and the dynamical part, entering the dynamical system of the model. 

It is easy to observe strong similarities of the model of the material point here 
with the geometrical formalism of homogeneous thermodynamics ( |161 126] . etc.) 
Exploiting these similarities, in the second part of this work we identify the phase 
space of a material point with the extended (by time) version of the thermodynam- 
ical phase space V of Caratheodory-Herman-Mrugala and the entropy form - with 
the dynamical part of the contact form in this space. This allows us to formulate 
the geometrical theory of material point and the corresponding entropy form in 
terms similar to those of the homogeneous thermodynamics, |26j . In particular, 
closeness condition of the entropy form is reformulated as the requirement that the 
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admissible processes curves belongs to the constitutive surface S of the model. Our 
principal result is the description of the constitutive surfaces of the material point 
model as the Legendre submanifolds Eg (equilibrium submanifolds of homogeneous 
thermodynamics) of the space V shifted by the flow of Reeh vector field. This shift 
is controlled, at the points of Legendre submanifold S5 by the entropy production 
function a. 

Ideally, the construction of a material point model (MP-model) in specific phys- 
ical circumstances should start by specifying the basic state space of physical fields 
and their spacial derivatives (gradients, divergences) whose time evolution one 
would like to study using the MP-model (notice that in that respect the MP-model 
is similar to the "Extended Thermodynamics", [31 ). Then one would like to deter- 
mine the dynamical system for these variables and, if necessary , to complement it 
with the constitutive relations closing the system. All this should be done in such 
a way that the energy balance law and the II law of thermodynamics expressed 
in the Clausius-Duhem inequality would be satisfied in a natural (for MP-model) 
form. In the Conclusion we mention some possible directions of this development, 
leaving its realization to the future work. 



Part I. Coleman-Owen model of thermoelastic point, entropy form and integrability. 

In this part we define and analyze the entropy form in the model of material point 
(MPM) suggested by B.Coleman and R.Owen, [7] and later on studied in numerous 
works [3 [9], etc. We will present our analysis on the example of thermoelastic 
material point used by B.Coleman and R.Owen as the basic model system. 

2. Entropy form of a thermoelastic system. 

In this section we remind the construction of the entropy form introduced 
in ([7]) and studied in [9]. We start with a balance of entropy of a continuum 
thermodynamical system in the form 

s + -V-Js = S, (2.1) 
P 

where s is specific entropy density, Js is the entropy flux and S is the entropy 
production that, due to the II law of thermodynamics, is nonnegative. Thus, we 
assume that the entropy supply is zero (the system is adiabatically isolated ). 

We also admit the relation between the entropy flux J5 and the heat flux q in 
the form 

3s = 0~\ + k, (2.2) 
where k is the extra entropy flux that will be taken to be zero in this section but 
will appear in more complex situations below. 

Plugging in the last expression, we rewrite (j2.ip in the form 

s + - V • (6l-lq) = s + -J- V • q + iq • VO'^ = E. (2.3) 
P Sp p 

Elastic material body is considered as a 3-dim material manifold Ad^ embedded 

at the time t into the physical (euclidian) space E^ by the diffeomorphism (pt : M ^ 

E^ . Deformation of M is characterized by the deformation gradient Fj — ^r)- 
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We admit the internal energy balance for a thermoelastic material point in the 
form 

pe = - V- q = cr:D-V-q. (2.4) 

where e is the internal energy per unit of volume, D = (Vv)'' is the symmetrized 
strain rate tensor (v being the velocity at the material point m), cr*-' is the Cauchy 
stress tensor and is the work power of the stress. Here and below the symbol : 
is used for the contraction of tensors. 

Expressing V • q from the energy balance and noticing also that D = (F~^F)'' is 
equal to the symmetrical part of the velocity gradient tensor L = F~^F we finally 
present the entropy balance in the form 

Op p 

Here we have used the angular momentum balance equation in the form = a. 

Now, we refer this relation to a material element, i.e. small enough volume of 
a material to associate with it definite values of the state variables participating 
in the energy and the entropy balance equations with the configuration, stress etc. 
see [7j. 

For such an element the Coleman-Owen model suggests a state space in a way 
that presumably guarantees the completeness of the dynamical system for the vari- 
ables in its state space (not yet chosen!) and describing for a given physical situation 
the response of the material element to the exterior influence. 

In such a case the evolution of the material properties is described by the col- 
lection of scalar, vectorial and tensor functions of time 

(p(t), F(i), 0{t), ci{t),e{t),a(t), {V 9-^){t) , etc.) , 

forming the process x the chosen material point. Variables listed above are 
related by some constitutive relations determined by the properties of the corre- 
sponding material media. Their time evolution has to be determined by a dynamical 
system. It is imperative for the closeness of the model to determine all such rela- 
tions and to use them for the reduction of the dynamical system to as simple one 
as possible. It is also important for the construction of dynamical models of the 
material point, see, for instance p!4l [T5] . 

The infinitesimal entropy production along a process x is given by the integral 
of 1-form 



ds - -edt + —{a- : F"^ : F)dt + -q ■ \/0^^dt = 
6 6p p 

^ds- -de + —a- : F^^ : dF + iq ■ VQ-^dt. (2.5) 
Op p 

Second expression in (2.6) can be considered as the exterior 1-form {entropy 
action form, see [9l fT2]) 

ri' ^-de- —a -.F-^ -.dF --q-VO-'^dt (2.6) 
9 dp p 
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for the process "f{t) considered as a time parameterized curve in the appropriate 
state space, i.e. the space of variables (e, F, t) and other variables cither independent 
in the proper state space or determined by the constitutive relations. 

The entropy action at the material point along a process 7(t) from the time 
T = to r = i is defined (postulated) as the integral 

As(7,)= / r/= f\*fj'- (2.7) 

Here 7*77' is the puUback of 1-form r]' to the interval of time. As defined, As(7t) 
is the functional of the curve j{t) with r G (0,t). So far nothing guarantees that 
the entropy action is defined in the space of processes 7 connecting two points in 
the state space. In particular it is unclear when the entropy function is defined as 
the function in the space of variables associating with a material point, up to an 
arbitrary constant (due to a choice of initial point). Assuming that the state space 
is simply connected this question rely on the property of the form 77 to be closed, in 
fact a drj — then locally (and in a simply connected space, globally)there exists 
the potential U 

r] = dU, AL/(7) = [/(7(end)) - U{-f{start)), (2.8) 

i.e. the function U of state variables (including time) defined up to an arbitrary 
constant. Notice that the potential U is defined by the entropy flux only and 
as a result coincides with the entropy (up to a constant) only when the entropy 
production in an admissible processes is zero. 

If the potential U and the entropy function s both exist as the functions in the 
space where processes 7 are studied, and if at the initial moment we normalize the 
potential U by the condition U{'j{0)) = 5(7(0)), the difference 

a = ,s(7(<)) - C/(7(t)) (2.9) 
is equal to the entropy production a during the process ^{t),t G [0,i]. 

Remark 1. Coleman and Owen postulated fulfillment of the II law of Thermody- 
namics in the inequality form that is weaker then the equality (|2.9p 

S{-i{end)) ~ S{-f{start)) ^ f 7*77'. (2.10) 

This property of upper semi-continuity of the entropy function S (seej^. Sec. 10) 
leads to the restriction on the space of direction of admissible processes. Namely, 
inequality (2.10) can be rewritten in the form 

t 

7*(d5-77') ^ 0. (2.11) 

If one has equality in this relation, both "/(t) and the inverse process "f{—t) 
are thermodynamically admissible. But if the inequality in ()2.11|) is strict, then the 
inverse process is prohibited. Since any segment of a thermodynamically admissible 
process has to be admissible, this conclusion can be localized, i.e. at each point 
m in the state space M, there is a cone Cm C Tm{M) of thermodynamically 
admissible directions. So, the II law requires the existence of a field of tangent 
cones in the state space defining directions of admissible processes. In continuum 
thermodynamics this is remedied by the "Amendment to the Second Law", see 
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[?, [27]. In the material point model this leads to some interesting geometrical 
consequences (see below). 

That is why in the papers [51 [HI [IHl [H] the conditions for closeness of this form in 
different situations were studied rather then more abstract notion of upper potential 
that would require much deeper thcrmodynamical analysis. 

2.1. State space and thermodynamical phase space. As the next step, we 
would like to note resemblance of the form lo = ds^rj' with the standard Gibbs form 
of the homogeneous thermodynamics (see |S] or the next section). Only the last 
term in (2.7) is qualitatively different form the terms in the contact Gibbs 1-form. 
To remedy this difference we suggest to use the heat vector field Ti introduced by 
M.Biot (Bi). The vector field H is defined uniquely, up to a constant in time vector 
field, by the condition 

n^q. (2.12) 

Using this vector field in the expression (2.7) for the 1-form r/', we transform it to 
the 1-form (non-degenerate by variables e,F,Ti.) 

rj=lde-^(T: F^^ : dF -(^9-^) : dn. (2.13) 
Op p 

At that point we introduce the configurational space of a material point 

B = {e,F,7i} (2.14) 
and the thermodynamical phase space (TPS) of the variables 

V = {s- = e, - F, =n;pi= 9-\p2 = j^o- : F-\ps = -^{^70-').} (2.15) 

The exterior 1-form 

tJ = ds-r]^ds- (ide - -^-cr : F'^ : dF - i(V6l-i) : dH) (2.16) 
Op p 

is typically contact (due to the functional independence of the variables (q^ ,pj) 
and, therefore, defines the contact structure in the space V. 



3. Integrability (closeness) conditions. 

Here we invoke the closeness condition of the entropy form both in (e, F, t) and 
{e,F,H) variables. In the Part II below it is shown that the closeness conditions 
of the forms 77, 77' are special cases of the integrability conditions for the contact 
structure in the extended thermodynamical phase space. 

For the beginning we revisit these conditions of integrability for the form 77', 
obtained in ?? and solve them. We remind that the basic phase space is {e, F, f3 ~ 
— j;^{d~^)} and the form ry' is 

T]' = dF + e-^de+{q-l3)dt. 

a 

Closeness conditions of rj' are 
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(3.1) 



dp{e-') = 0, 

The first three equations have the general local solution of the form 

(e-^ = dM + ci{t), 

for a function [/(e,F, i) and arbitrary functions ci{t), C2{t). 

Using these expression in the second triplet of equations one sees that these 
three equations in p.ip are equivalent to the following presentation of the function 
(q-/3): 



dU 

(q-/3) = — +4(i)-F + c;e + C3(i), 



(3.2) 



with an arbitrary function c^it). 

The function U{e,F,t) is the time-dependent entropy form potential defined in 
the space of basic variables. These formulas determine the following constitutive 
relations 



-1 = a,c/ + ci(t), 

. -p-l _ dFU+C2{t) 

a,(7+ci(t))' 

_(q • /3) = dtU + 4(t) • F + c[e + csit). 



(3.3) 



Repeating the same arguments for the representation ry of the entropy form in 
the space of variables {e, F, Ti} we get the following 

Proposition 1. For the model presented above with the space of basic fields {e, F, Ti.} 
condition of integrability of the entropy form rj given by i2.13\) is equivalent to the 
existence of a potential U{e, F, 7t) such that the following constitutive relations hold 



(T : F-i = -p{d,U)dFU, 
V(6i-i) = -pdnU. 



(3.4) 



4. Partial integrability, dynamical equations and constitutive 

relations. 

Here we consider the dynamical system for the basic variables e,F,7t reformu- 
lating the dynamical system introduced by Coleman and Owen: 




(4.1) 
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Here 7 is just the notation for the derivative of /3 = — p ^. 
As the first step we replace third equation by the relation defining Ti.: 

F = LF, 

-p-imL-p-iV-q, (4.2) 

n = q. 

Before considering the full system of integrability conditions (|3.ip . let us look at 
the meaning of partial integrability of the entropy form; first for individual terms 
(which, of course, is always possible) and then for couples of these terms. 

More specifically, assuming the Fourier relation between the heat flux and the 
temperature gradient we have the relations 




Thus, if V6'"^ can be determined as a function of basic variables , one of the 
equations of the dynamical system will be closed. Looking at the form ry, see (j2.13p 
(or at the equation ([3?l])3 we notice the integrability of the third term, i.e. the 
possibility to write 

--(V6I-I) : dn = dn^dn 
P 

with a function ^(e,F,7^) that guarantees fulfillment of {V9^^) = pd-j-cS,, and 
allows to rewrite third equation in the system (4.2) in the closed form 

n = k-'p-^{e,F,n). (4.3) 

Thus, integrability of the third term in (|2.13p allows to close the third 
dynamical equation of the system (4.2) in the space of basic variables. 

Integrability of the first term of the form 77 is equivalent to the statement that 

for some differentiable function W{e,F,'H). 

Reversing this relation we get the constitutive relation in the form e = e{d, F, Ti.). 

Looking at the systems (|4.1I) and (|4.2p we see that the integrability of the first 
term in (|2.13p delivers the basic constitutive relation presenting internal 
energy e as the function of basic variables F,7t and temperature 9. 

Integrability of the second term of the form rj is equivalent to statement that for 
some function Z(e,F,?<), 

— a = dpZ : F. 

Op 

It expresses the Cauchy stress as the function of basic variables. As a result, 
integrability of the second term in (|2.13p allows to determine the elas- 
tic part of the evolution of the internal energy e in the space of basic 
variables. 

If we allow the form 77 to be integrable in the first two terms, functions W, Z 
should coincide up to an arbitrary additive term, i.e. A(F): Z = W + X{F) and, 
using the expression (j4]) for 6^^, we get 
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-a- : = (d^Wr^drW. 
P 

We will get similar simplifications assuming integration other couples of terms in 
the expression (|2.13p for the form -q. Finally, if we assume the full integrability of 
77, we get the system (4.2), in terms of a potential J7(F, e,?t) of the entropy form 

= LF, 

=.-^|^:L-p-iV.q, (4.4) 



H — k ^ pWj 



'e dF 
PL 

on- 



Additionally we get the constitutive definition of temperature through the po- 
tential U as: 



Canonical equation for F 



F = LF 

has pure kinematical meaning. 

Thus, to close the dynamical system the for variables (e, F, Ti.) one has to deter- 
mine, by some constitutive relations or phenomenologically, the fields L = Vv and 
V • q. Another way would be to include them in the list of dynamical fields and 
construct additional dynamical equations for them in the spirit of rational extended 
thermodynamics, see (31]. The first way led Noll to his definition of generalized 
processes [3D]. The arguments presented above show that the constitutive relations 
required in the definition of Noll reduces to the heat propagation constitutive rela- 
tion (Fourier, Cattaneo, etc.) 



5. Entropy form in deformable ferroelectric crystal media 

Let us now recall the phenomenological model of an elastic deformable ferroelec- 
tric crystal medium worked out in [lI]-[13ji[Mj in a suitable Galilean quasi-static 
approximation. We assume that the medium is formed by n molecular species, each 
one of them giving rise to a field of electric dipoles. The total polarization per unit 
of mass is given by: 

- - (5.1) 



P 

where Ppoi is the total polarization (per unit of volume) 

The vector tt is assumed to satisfy the following balance equation 

Iii = £ + ^E + p-\\7 -^E), (5.2) 

where / 7^ is an " inertia constant" (which in the following will not be restrictive 
to let be equal to one), £ is an external electric field calculated in a comoving frame, 
■^E is a vector field, called local electric field, accounting for the interaction between 
the polarization of different species with the crystal lattice, V E is the divergence 
of a rank-two tensor ^E, called local electric field tensor, which accounts for the so- 
called shell-shell interaction and finally ^E is responsible for the typical ferroelectric 
ordering. The master equation (|5.2p fixes the behavior of the state variable tt in 
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terms of the electric fields ■^E and ^E. Fields and play the role of internal 
variables. This equation resembles Newton equation law of motion. 
The internal energy balance has the form: 

pe-P(,)-V-(q-7') (5.3) 
where 'P is the Poynting vector and the work power of internal forces is given 

by 

= cr : L - p^E : TT + : Vtt (5.4) 
By using the relation L = FF^^ we can rewrite (j5.4[) as: 

= cr : F"i : F - p^E : iV + -^E : VTi-. (5.5) 

Then, the dynamical system describing evolution of the state fields of the system 
[l2] is the following: 

F = LF 

e = (T : F-i : F - p^E : tV + ^E : - p-^V • (q - P) 
H = q 

< tV = u (5.6) 
u = £: + -^E + p-i(V-'^E) 
Vtt = Vu 

,Vu = V • Jvu + CTVu, 



where Jvu and ctvu are both phenomenological quantities representing, respec- 
tively, the current and the source terms associated with Vu. The introduction of 
the variable u = tt has been used so to obtain a first order dynamical system. 

We finally present the entropy balance ()2.3p in the form: 

s-e^^k+pB~^a : F"^ : F-6I"^^E : 7T + (p)"^^E : V7T-p~^q-V6l~^-p"^V-(6l~^7'+k) = H 

(5.7) 

where the general relation (j2.ip for the entropy flux has been assumed. 
Then, the infinitesimal entropy production along a process is given by the integral 
of the following 1-form: 

do = ds - e^^de + pO^^a : F"^ : dF - e^^^E : (ivr + 7r+ 

+ p^i^E : Vd - p-^ci ■ Ve-^dt - p-^V ■ {Q-^V + k)dt (5.8) 

In such a case the evolution of the material properties is described by the fol- 
lowing functions of time forming the state space 

B = (e,F,7Y;7r, V7r,7T-, V7r;t). (5.9) 

The vector field has been introduced above. 

The right side of (|5.8p can be considered as the exterior 1-form: 
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Tj = e-^de - pQ-^cr : F"^ : dF + 6'-^^E : dn - pO-'^^E : dVn + (3 ■ d'H+ 

- p-'^{9-^\7 -V + V ■k)dt, (5.10) 

remembering the definition f3 = . 

We notice that the Poynting vector together with the k quantity determine 
outside contributions in the entropy production. These arguments lead us to the 
introduction of the large space of the variables: 



B 



sUBU 



-1 



^E,-V6'"\-V-k,--V-7'] 



(5.11) 



The exterior 1-form u = ds — rj defines the contact structure in the space B. We 
notice that the differential of the form u is given by 



n = dLj^ -d-q. (5.12) 

If r] is closed then, locally, r] — dU for some functions U £ C°°{B). Assuming 
that the state space B is simply connected the question of the existence of the 
entropy as function of the state variables rely on the property of the form 77 to be 
closed, i.e. drj = 0. This gives the reason to analyze the conditions for closeness of 
this 1-form. 

Closeness of the form gives the potential ?7(e, F, tt, Vtt, ?t, t) on B so that 

-(p6l)-icr-F-^ = 5fC/, 

(5.13) 

-{pey^^E^d^^u, 

J3\7 - V - p-^V ■k = dtU. 

By using the identity ^-^V • P = V • - boldsymbolV ■ WO''^ the last 

equation of the above system takes the form 



T ■ 13- p-^V -{e-^r + 'k) =dtU. (5.14) 

Finally, if we assume the full integrability of 77, we get the system ()5.6p in terms 
of potential U of the entropy form as: 
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F = LF, 

e = edpU : LF - pOd^U ■ u - pOd^^ : Vu - p-^V ■ q - 6l(pV • k + dtU), 
n = k-^pdnU, 

< TT = U, 

u = 5 + 0d^U - p-\V ■ epd^^U), 
V'tt = Vu, 

Vu ^ V • Jvu + crvu- 

(5.15) 

In order to close this system of equations in the space of basic fields B one would 
have to use constitutive relations for the following entries in the system: 

(1) The terms L, V • q as in the simple model [7], 

(2) the term V-k which would be present if mixed dissipative processes would go 
in the system, that one usually is determined using a dissipative potential, 

(3) gradients of variables appeared in V • 9pd\jTrU), 

(4) flux and production terms in the last equation. 

It is hardly possible to do this in some regular and relatively simple way. This 
problem is similar to the problem of closeness of systems of equations for momenta 
in Statistical Mechanics or Rational Extended Thermodynamics. 

It seems more reasonable to lift the problem to the larger space, i.e. the ex- 
tended thermodynamical phase space B where several models of thermodynamical 
processes were suggested, [HI [H [15] , etc. 

To prepare the framework for such development, we present, in the next part of 
this work, the geometrical scheme introducing the entropy form, its integrability 
(closeness) into the conventional thermodynamical phase space, extended, whenever 
necessary, by the time t added to the list of extensive variables. 



Part II. Contact geometry of entropy form. 

6. Contact structure of homogeneous thermodynamics. 

In this section we briefly recall the standard contact structure of homoge- 
neous thermodynamics in the thermodynamical phase space introduced by C.Caratheodory 
and developed by R.Hermmann and R. Mrugala f [51 fTB]'). 

A phase space of the homogeneous thermodynamics (thermodynamical phase 
space, or TPS) is the (2n-|-l)-dimensional vector space V = endowed with 

the standard contact structure ([21 [16]). Contact structure is defined by the 
(contact) 1-form such that exterior product of "d and n copies of its differential 
dd is nonzero (2n -|- l)-form: 

§ ^d^ ^...^d^ 

By D'Arbois Theorem, [l], there is a choice of coordinates (local in a general man- 
ifold and global for the standard contact structure) (z; (g^,pi), . . . (g"',p„)) such 
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that the takes the form 

n 

d^dz^^Pidq\ (6.1) 

i=l 

The horizontal distribution D = Ker{d) of this structure is generated by two 
famihes of vector fields 

D =< (9p,, ^q^ +pids > . 

The 2-form 

n 

n = d-d = - ^dp, A dq' 

i=l 

is a nondegenerate, symplcctic form on the distribution D. 

The Reeb vector field, uniquely defined as the generator C of the 1-dim charac- 
teristic distribution ker{dd) satisfying 'i?(C) — 1, is simply 



7. GiBBS SPACE. LeGENDRE SURFACES OF EQUILIBRIUM. 

Concrete thermodynamical systems are determined by their constitutive rela- 
tion, which, in their conventional form determine the value of a thermodynami- 
cal potential z = E{q^) as the function of n (extensive) variables of a D'Arbois 
canonical coordinate system (z; {q^,pi), . . . {q"',Pn))- Dual, intensive, variables are 
determined then as the partial derivatives of the thermodynamical potential by the 
extensive variables: pi = 

Geometrically, a constitutive relation is determined as a Legendre subman- 
ifold (maximal integral submanifold) E^; of the contact form i?. Locally a Legendre 
submanifold S is determined by a choice of canonical coordinates (z; {q^ ,pi), . . . ((f' ,p.. 
such that an open subset C/ C E projects diffeomorphically to the space X of vari- 
ables q^ . In terms of these coordinates E is defined in the open domain U as follows: 

dE 

E = {{z, q,p) e P\z = E{q'),p, = (7.1) 
more about local presentation of Legendre submanifolds and their properties see 

mm- 

Space G of variables z, g', i = 1, . . . ,n is, sometimes, named the Gibbs space 
(bundle) of the thermodynamical potential E{q^). Thermodynamical phase space 
(■p,??) (or, more precisely, its open subset) identifies with the first jet space 
J^iY X) of the (trivial) line bundle tt : G X. Projection of Eg to the Gibbs 
space Y is the graph F^; of the fundamental constitutive law E — E{q^). 

Another choice of the thermodynamical potential together with the n-tuple of 
extensive variables leads to another representation of an open subset of TPS V as 
the 1-jet bundle of the corresponding Gibbs space. 

The most commonly used thermodynamical potentials are: internal energy, en- 
tropy, free energy of Helmholtz, enthalpy and the free Gibbs energy. 

On the intersection of the domains of these representations, corresponding points 
are related by the contact transformations (see jllll])- 

Example 1. As an example of such a thermodynamical system, consider the van 
der Waals gas - a system with two thermodynamical degrees of freedom. Space V 
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is 5-dimensional (for 1 mole of gas) with the canonical variables {U, (T, S), {—p, V)) 
(internal energy,temperature, entropy, -pressure, volume), the contact form 

i3 = dU - TdS + pdV, 

and the fundamental constitutive law 

U{S, V) = {V-b)^e^ -y, 

where R is the Ridberg constant, cy is the heat capacity at constant volume, a, b 
are parameters of the gas reflecting the interaction between molecules and the part 
of volume occupied by molecules respectively, see [1] . 

8. Extended thermodynamical phase space, its contact structure. 

Introduce the extended thermodynamical phase space T^^m+i j^jrj'pps) 
This space contains m physical fields = 1, . . . ,to that may include, together 
with basic fields (temperature, density, polarization vector, etc.) also their space 
gradients. As an element new in comparison to the usual thermodynamical phase 
space in V, time t may be included as the m-st variable = t. In addition, 
m variables pi dual to (including, possibly, p,„ dual to the time variable = t) 
and the thermodynamical potential (entropy in our work) s are considered as the 
variables in the space V. It is convenient to consider both types of situations - 
where time is included as an independent variable and where it is not (see Section 
2 where the field Ti was introduced). 

Remark 2. We change notation of the coordinate z of the canonical coordinate 
system into s because in the case of MT-model, we will consider entropy s as the 
thermodynamical potential. 

ETPS V will be endowed with the standard contact structure 

m 

i}^ds-J2p^dq' (8.1) 

i=l 

As above we denote hy D — Ker{-d) the horizontal distribution of the contact 
structure (P,"!?) and by ^ the corresponding Reeb vector field. 

Definition 1. (1) An extended constitutive surface (ECS) is the m-dim 
submanifold E C "P such that the restriction of the contact form -i? to S is 
exact: 

i9|s = rf(s-f/) (8.2) 

for a function U E C°°(I]). 
(2) If S C P is an extended constitutive submanifold and U - function in the 
definition above, the function 

a — s — U, 

defined on the surface E is called the entropy production potential. 

Remark 3. It is convenient to write the function in the right side of (|8.2p as s — J7 
because this U coincide with the potential of the entropy form rj from the Part I. 
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Proposition 2. For an arbitrary extended constitutive surface E there exists a 
Legendre submanifold So of the contact structure d cmd, a, function a S C°°(S]o) 
such that the ECS S is obtained from Legendre submanifold Sq by the shift by the 
flow of the Reeb vector field (: A point {s,q^,Pi) G So is shifted with the value of 
parameter equal to (—a{s,pi,q^): 

Y. = {exp(-cr{s,p,q)C)(m), m = (s,p, g) £ Sq.} (8.3) 

Proof. If (S, U) are respectively an ECS and the corresponding potential, extend 
the function U out of the surface S to some smooth function in V. The form 
— Y^^QPidq'' in V vanishes on the submanifold S. 

Consider the mapping (f> : {s,pi,q^) — *■ {s + a,pi, q'') in a neighborhood of surface 
S, Let (f){T,) be the image of surface S under this mapping. Notice that the puUback 
of contact form ^? under the mapping ^ is: 

m 

(l>*^ = d{s-U)-^Pidq\ 
Let S, G Tx{T,) be a tangent vector to the surface S at a point x. Then 

m 

=< d{s -U)- J2Pi(^<l'^ ^ >=< ^ >=< ^' -^^^ > • 

Therefore, surface Sq = </>(£) is Legendre surface of contact structure {P,'d)- 
Therefore, the ECS surface S is obtained by the deformation of a Legendre sub- 
manifold So by the flow mapping exp{—(jQ of the Reeb vector field C, at the value 
of parameter —a. □ 

Remark 4. Function U is defined up to addition of a constant. Another choice 
of U leads to the shift of the ECS by this constant in the direction of variable s, 
i.e. along the trajectory of the Reeb vector field C,. Correspondingly, the Legendre 
submanifold Sq is shifted. As a result, we came up with the (rn+ 1) — dimensional 
submanifold A foliated by the shifts of Legendre submanifold Sq ECS S[/ and, 
transversally, by the phase curves of the Reeb vector field. 

Corollary 1. With any admissible process % : T — > S (see below the definition) 
there is related an uniquely defined reversible process xo '■ T ^ T,o defined by the 
condition: 

Xo{t) = expiaix{tmox{t), (8.4) 

where 

a^s-U (8.5) 
is the entropy production potential. 

Denote by tte : S — > X the projection of ECS to the base X and by tteq - 
corresponding projection for associated equilibrium surface Sq. 

If this projection is invertiblc on some open subset W C T,, denote by : 
U ^ W C corresponding inverse mapping. We will call mapping js - the 
characteristic embedding of S. 

Using the embedding j one can directly relate the closeness of entropy form rj 
studied in Part I with the requirement of integrability of the contact form along 
the admissible dynamical processes x{t)- 
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Proposition 3. Entropy form 77 is closed on the base space X if and only if 

11 = f§ {= dU) (8.6) 

where j : X ~> V has the extended constitutive surface as its image. 

A conventional way to specify such an extended constitutive surface is to identify 
the ETPS (P, I?) with the 1-jet space of the Gibbs bundle - a hne bundle G ^ X 
over the space X of variables {q'^,t),see [16], (or simply in the time-independent 
representation): 

G = RsXX (8.7) 



X 

A choice of a section of the bundle tt - an entropy function s = S{q'^,t) allows 
to form the Legendre submanifold S5 — j^(S)(X) - image of the space X under 
the 1-jet section j^{S) of the bundle V G ^ X . Legendre submanifold Ss 
constructed in such a way projects diffeomorphically to the space X. 

Choose next a function cto G C°°{X). Lift this function to the surface S5 to get 
the function 

cr = 0-5 o ttiIes- (8.8) 
Now we define the extended constitutive surface 

Ss,a - {z* - expi-a{z))z\z £ Eg}. (8.9) 

In this case the mapping j{q) = exp{—a{q)) o j^(S){q) is the characteristic em- 
bedding for ECS Ss,^- 

It follows from the local description of Legendre submanifolds of contact struc- 
ture ([rif^ that the following statement is valid 

Proposition 4. Any extended constitutive submanifold T, G V with the set of basic 
variables locally has the form Ss.ct for two functions S*, cr G C°°{X). 

Remark 5. For a general Legendre manifold, where the function ("potential") 
defining the variable of z depends on q^, i G I',Pj,j G J for some decomposition 
[1,to] = I U J of the set of indices from 1 to m, one can modify this definition 
accordingly to define corresponded shifted Legendre submanifold. 

Notice, that since a general Legendre submanifold of the contact manifold {V, 9) 
might have singularities or be given by multivalued function 5", a general extended 
constitutive surface S and its projection to X may have singularities (Legendre 
singularities, see ([!]). A standard example of this kind is the projection of the 
constitutive surface of the van der Waals gas to the pT plane, see [4]. 

Let X : T — > X be a curve defining the process (evolution of state). Combining 
it with the constitutive mapping js^„ of a ECS T,s,a we get the curve x the space 

r. 
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Definition 2. A t her mo dynamically (TD) admissible process (with the en- 
tropy function S and the entropy production potential a) is the curve x '■ ^ ^ 
such that for its lift x : T — > T,s,a given by 

x{t) = jsAx{t)), 

one has 

imit)), X'{t)) = dd(x'(t)) - {da(x(t)), x'it)) ^ 0, for all t. (8.10) 

Thus, the calculation of the entropy production in this model can be done directly 
on the space X of basic fields. Change of entropy along the process x during the 
time interval {tQ,t\) is equal to 

As|*; = AC/ + Aa = (C/(ti) - U{ta)) + {a{h) - a{to)). (8.11) 

Resume: In a case of the closed entropy form r], entropy form model is (at 
least locally) presented by two potentials - {U, a) or {S = U + a, a) where U is the 
potential of the form r], S is the entropy density and a is the entropy production 
potential. Potential U (or S) is defined up to addition of a constant. 

Dynamical evolution is presented by a thermodynamically admissible curves in 
the space X of basic fields or by a the TD-admissible curve on the extended con- 
stitutive surface T,s,a in the extended thermodynamical phase space {V, i?). 

9. Entropy form as a flat connection in the Gibbs bundle. 

Gibbs space G is endowed with (local) coordinates {s;q^,i = 0, . . . , m) and the 
corresponding frame dt, dqi,ds together with the corresponding coframe dg', ds. 

Having the constitutional surface 'Es.a and the corresponding Logcndrc subman- 
ifold Sg available we may extend to the (m-l-l)-dim submanifold As, possibly 
with singularities (or locally, without singularities) by applying to the points of S5 
the flow of the Reeb vector filed C: 

As = {exp{{s - S{q))Ojsms,q) e G}. (9.1) 

This definition of A5 allows to define the smooth mapping 

A : G ^ As : X{s, q) = exp{{s - S{q))C)js{q), 
which, away from the possible singularities of Lagrange submanifold Ss is correctly 
defined. 

In this situation we can associate with the contact 1-form 9 = ds — ^ ■ Pidq^ in 
the ETPS V the 1-form 

m 

u) = X*e = ds-r] = ds-^Pi{q)dq\ (9.2) 

i=0 

Constructed 1-form defines the projector 

n,(0 = {ds - m)ds, C G T{G) (9.3) 
to the vertical subbundle V{G) C T{G) of the tangent bundle of space G: 

This projector defines the connection iv on the Gibbs line bundle. 
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Horizontal space of this connection is defined by the condition n„(^) = 0. It 
consists of the tangent vectors ^ = £,^dt + C^qi + ^''ds for which 

so, that 

Hor{cu) = {< dt+Pods,dg^ + pA >}. (9.4) 
Coefficients of the projector form do not depend on s, thus the connection defined 

in this way is the Hnear connection. We wiU call it ?7s since it is defined by the 

constitutive relation E. 

Curvature of the connection 775] is zero 

ft = D„Tjs = -drjs{nh-,Ilh-)) = (9.5) 

since rj^, = dS. 

9.1. General connection in Gibbs bundle as an entropy form. Consider a 
general connection in the Gibbs line bundle for a situation where no t is present in 
the list of canonical variables (or simply take t — q™~^^, i.e. defined by the 1-form 

V ^ Pi{s,q)dq\ (9.6) 

where functions pi{s,q^) are arbitrary. Connection form defining projection to the 
vertical subbundle (ds) C T{G) is 

Lu = ds — rj. (9.7) 

Horizontal vector fields have the same form (|9.4p as above but with ^0 = 0. 
The curvature of this connection is determined by the curvature form f2 = Di^u — 
duj{Ph- ,Ph-) m] Since 

dri = Pi.sds A dq^ + Pijdq^ A dq\ (9.8) 
then on the couples of basic horizontal vector fields dgi + Pidg, 

il{dqi + Pids, dqj + p-jds) = {P],sPi - Pi.sP]) + (Pj,g- " Pi.qi)- (9.9) 

In this general case integrability conditions for the form 77 are not fulfilled and 
entropy potential U is not defined even if the entropy form is independent on 
variable s. 

Remark 6. It would be interesting to see, possibly, on examples, the meaning of 
curvature in the case where entropy form is only partly integrable. 

10. Entropy form geometrically: conditions of integrability. 

In a number of works, see [H [9l [lOl [12] and the literature cited thereon, the close- 
ness conditions of the entropy form were calculated for different thermodynamical 
systems. Here we show that these conditions take simple abstract form in terms 
of the contact structure of ETPS V. In these form the integrability conditions 
may be useful for writing down and studying similar conditions for specific physical 
systems. 

Consider the 1-form in the extended base space with variables (g* , t) 

r/ ^p,{q\t)dq' -H{q\t)dt. (10.1) 
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Let the form 77' be closed: drj' = 0. Then 

= {p^^qjdq^+p^fidt)Adq''-{H^gJdq^+H^tdt)Adt ^ -{pi^^j -pj^gi)dq^ Adq^ -{p^^o+H^qi)dq^ Adt. 

Thus, closeness condition for the 1-form 77' is equivalent to the fulfillment 
of the following relations: 

[Pi^t + H^q, = 0, Vi. 
and, locally, is equivalent to the existence of the potential U: 

where the function U (g*, t) is defined uniquely up to adding of an arbitrary constant. 

Let s {q{s) , p{s)t{s) be a parameterized curve in the phase space. Then, we 
have the relation between rates of change of p and q variables: 

dsPi = U,q.qjdsq^ + U^g^t^st■ 

In particular, for a curve t {t, q{t),p{t)) we have 

dtP^^U^q.q^dtq' +U^q^t. (10.4) 

These relations are of geometrical nature and has to be compatible with any dy- 
namical evolutional model of the material point. 

Remark 7. Hessian matrix Rij = U_qiqj of the entropy potential U defines the 
thermodynamical metric.^6\ , coinciding with the Ruppeiner metric [29j (defined 
by the Hessian of the entropy S) on the constitutive surface Pi = §^ where the 
entropy production a is constant. Equation (10.4) represents the relation between 
tangent and cotangent vectors to the process curve on the base X defined by this 
metric. 

Condition of non-degeneracy of this metric is: 

det{Uq.q,)^deti^)^0, 

i.e. represents local condition of invertibility of the godograph transformation — ^ 
Pi- 



ll. Conclusion. 

In this work we analyzed integrability (closeness) conditions of the entropy form, 
introduced by Coleman and Owen in the model of material point. 

Considering first the model of thermoelastic point and then the thermoelastic 
dielectric material, we showed that the entropy potential that is defined for the 
closed entropy form determines some terms in the assumed dynamical system of 
corresponding model and produces the constitutive relations for temperature and 
its gradient. 

On the other hand, dynamical systems for these systems stayed unclosed and 
the arguments are presented that in the space of basic fields it is impractical if not 
impossible to close this system. 
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That is why as an alternative to the basic space of fields, we suggest to use 
extended thermodynamical phase space (ETPS) ("P, t9) with the canonical contact 
structure "d similar to one studied in homogeneous thermodynamics. 

In the second part of work we prove that the integrability condition of the entropy 
form is equivalent to the conditions that processes are confined to the "extended 
constitutive surface" S of ETPS. 

Structure of such surfaces is determined in terms of Legendre submanifolds Sq of 
the standard contact structure of ETPS P and additional function S having sense of 
the entropy production potential. As a result, a material point model in the ETPS 
with the integrable entropy form is defined (at least locally) by two functions - 
entropy S - (function of basic fields, defining the Legendre submanifold Sg) and 
the entropy production potential a on such a submanifold defining the shift along 
the phase curves the Reeb vector field that produces the extended constitutive 
surface S. 

Thermodynamically admissible processes are defined as the time parameterized 
curves on the surface T, along which the value of entropy production potential is 
increasing. This geometrical model of the material point represents, in our opin- 
ion, an adequate framework for construction of a thermodynamically admissible 
dynamical model. 

Next step in the development of a material point model(s) would be to study the 
compatibility of known geometric models of irreversible thermodynamical processes 
- Lagrangian systems with dissipative potential, [20 j . metriplectic systems, 
and the gradient relaxation processes of H. Haslash,[T4l [15] with the geometrical 
material point model. This problem will be studied elsewhere. 

References 

[1] V.Arnold, Mathematical Methods of Classical Mechanics, Springer, 1989. 

[2] V.Arnold, A.Givental, Symplectic Geometry, in "Dynamical Systems IV", Springer, 1988. 

[3] M.Biot, Variational Principles in the heat propagation, 1961. 

[4] H.Callen, Thermodynamics. 

[5] CCavatheodory, Unterschungen uber die Grundlagen der Thermodynamyk, Mathematishe 
Annalen, 67, 1909, 355-386. 

[6] V.Ciancio, M.Dolfin, M.Francaviglia, P.Rogolino, Thermodynamic Transformations in Mag- 
netically Polarizable Undeformable Media,J.Non-Equilib. Thermodyn.,2001, v. 26, pp. 255- 
267. 

[7] Coleman, B.D., Owen,R.D. Mathematical foundations for Thermodynamics R.D.,Arch. Rat. 
Mech. Anal.,54,pp.l-104, 1974. 

[8] M.Dolfin, M.Francaviglia, P.Rogolino,/! Geometric Perspective on irreversible Thermody- 
namics with Internal Variables, J.Non-Equilib. Thermodyn., 1998, v. 23, pp. 250-263. 

[9] M.Dolfin, M.Francaviglia, P.Rogolino, A geometric model for the thermodynamics of simple 
materials. Periodica Politechnica Ser. Mech., Engr., v. 43, n. 1,29-36, 1999. 
[10] M.Dolfin, M.Francaviglia, L.Restuccia, Thermodynamics of Deforamable Dielectrics with a 
Non-euclidian Structure as Internal Variable, Technische Mechanik, B.24, H.2 (2004), 137- 
145. 

[11] D.Eberard, B.Maschko, A. J. van der Shaft, An extension of Hamiltonian systems to the ther- 
modynamical phase space: Toward a geometry of nonreversible processes. Rep. Math. Phys., 
V.60, no.2, 2007, pp. 175-198. 

[12] M.Francaviglia, L.Restuccia, P.Rogolino, Entropy production in polarizable bodies with in- 
ternal variables, J.Non-Equilib. Thermodyn., 2004, v. 29, pp. 221-235. 

[13] Godbillon, C, Geometric Differentielle et Mechanique Analytique, Hermann, Paris, 1969. 

[14] H.W.Haslash, Geometric Structure of the Non-eguilibrium thermodynamics of homogeneous 
systems. Rep. Math. Phys., v.39, no.2, 1997, pp. 147-162, 



GEOMETRY OF THE ENTROPY FORM. 



21 



[15] H.W.Haslash, A non- equilibrium therrnodynamical geometric structure for thermoviscoplas- 
ticity with maximum dissiption, Int. J. of Plasticity, 18 (2002), 127-153. 
R. Hermann, Geometry, Physics and Systems, Dekker, N.Y., 1973. 

S.Kobayashi, K.Nomizu, Foundations of Differential Geometry, Wiley (Interscience) , N.Y., 
Vol.1, 1963; Vol.11, 1969. 

J.Marsden, T.Hughes, Mathematical Foundations of Elasticity, Dover, N.Y.,1983 
G.Maugin, Material Inhomogeneities in Elasticity, Chapman and Hall, 1993. 
Maugin, G., Internal Variables and Dissipative Structures, J. of Non- Equilibrium Thermo- 
dynamics, 5 (1990), 173-192. 

Maugin G.A., Deformable Dielectrics- 1- Field Equations for a Dielectric Made of several 
molecular species, Arch. Mcch. (Poland), 28,(1976), 676-682. 

Maugin G.A., Deformable Dielectrics-II- Voigt's intramolecular force balance in elastic di- 
electrics, Arch. Mech. (Poland), 29,(1977), 143-159. 

Maugin G.A., Deformable Dielectrics-III- A model of interactions, Arch. Mech. (Poland), 
29,(1977), 251-258. 

Maugin G.A., Pouget J. Electroacoustic equations for one-domain ferroelectric bodies, 
J.Acoust.Soc.Am, 67,(1980), 588-601. 

, Mrugala, R, Nulton, J.,Schon, J., Salamon, P., Contact Structure in Thermodynamic The- 
ory, Reports on Math. Phys., v.29, No.l, pp. 109-121, 1991. 

R. Mrugala, Geometrical Methods in Thermodynamics, in Thermodynamics of Energy Con- 
version and Transport ed. S.Sieniutycz, A.de Vos., Springer, 2000, pp. 257-285. 
Muschik,W., Aspects of Non- Equilibrium Thermodynamics, World Scientific, Singapur, 1990. 
S.Pickarski, GALILEAN-INVARIANT FORMULATION OF THE FLUID MEGHANICS, 
IFTR REPORTS 7/2007. 

G.Ruppeiner,Riemannian geometry in thermodynamic fluctuation theory, Reviews of Modern 
Physics, V.67, n.3, 1995, pp.605-659. 

C. Truesdell, W. Noll, The Non- Linear Field Theories of Mechanics, 2nd ed.. Springer, 1992. 
l.MuUer, T.Ruggeri, Rational Extended Thermodynamics, 2nd ed.. Springer Tracts in Natural 
Philosophy, Vol.37, 1998. 

Department of Mathematics, University of Messina, Messina, It. 
E-mail address: dolf ln{Sdlpmat .unlme . it 

Department of Mathematics and Statistics, Portland State University, Portland, 
OR, USA 

E-mail address: sergeOmth.pdx . edu 



Department of Mathematics, University of Messina, Messina, It. 
E-mail address: lrestSdipmat.unime.it 



